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^^ I In this article, we give a polynomial algorithm to decide whether a given permutation a is 

^P ■ sortable with two stacks in series. This is indeed a longstanding open problem which was first 

introduced by Knuth in [I]. He introduced the stack sorting problem as well as permutation 
patterns which arises naturally when characterizing permutations that can be sorted with 
one stack. When several stacks in series are considered, few results are known. There are 
>-^ ■ two main different problems. The first one is the complexity of deciding if a permutation 

T^ ', is sortable or not, the second one being the characterization and the enumeration of those 

C^ '_ sortable permutations. We hereby prove that the first problem lies in P by giving a polynomial 

algorithm to solve it. This article strongly relies on [3] in which 2-stack pushall sorting is 
defined and studied. 

O ■ 1 Notations and definitions 

\o 

5^ ' Let / be a set of integers. A permutation of / is a bijection from / onto /. We write 

• . a permutation cr of / as the word a = a\a2---On where Oi = o"(ii) with / = {ii...in} 

(<— ^ I and ii < 12 < ■ ■ ■ < in- The size of the permutation is the integer n and if not precised, 

cn ' I = [1 . . . n]. Notice that given the word (Ti(T2 ■ ■ ■ ctu we can deduce the set / and the map 

a. For any subset J of /, a^j denotes the permutation obtained by restricting a to J. In 

particular the word corresponding to a\j is a subword of the word corresponding to a. 
k> ' Let's recall the problem of sorting with two stacks in series. Given two stacks H and V 

^ ■ in series -as shown in Figure [1]- and a permutation a, we want to sort elements of a using 

the stacks. We write a as the word a = (Ji(J2 ■ ■ -CTn with Uj = a[i) and take it as input. Then 

we have three different operations: 



C^ 



p which consist in pushing the next element of a on the top of H. 

A which transfer the topmost element of H on the top of V . 

jjL which pop the topmost element of V and write it in the output. 
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Figure 1: Sorting witli two stacks in serie 



If there exists a sequence w = wi . . .Wk oi operations p, A, /i that leads to the identity in 
the output, we say that the permutation a is 2-stack sortable. In that case, we define the 
sorting word associated to this sorting process as the word w on the alphabet {/?, A, //}. Notice 
that necessarily w has n times each letter p, A and /i and k = 3n. For example, permutation 
2431 is sortable using the following process. 
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This sorting is encoded by the word w = ppXpXpX^X^fifj,. We can also decorate the word 
to specify the element on which the operation is performed. The decorated word for w and 
2431 is vj = P2PiXiP2X2,piXipiX2p.2P^Pi- Note that we have the same information between 
{(J^w) and w. Nevertheless, in a decorated word appears only once each letter pi^Xi or pi. 
The decorated word associated to (ct, u)) is denoted w" . 

Of course not all permutations are sortable. The smallest non-sortable ones are of size 7, 
for instance a = 2435761. 

When only one stack is considered, there exists a natural algorithm to decide whether a 
permutation is sortable or not. Indeed, there is a unique way to sort a permutation using 
only one stack, and a greedy algorithm gives a decision procedure. For two stacks in series, 
a permutation can be sorted in numerous ways. Take for example permutation 4321. Each 
element can be pushed in either stacks H ox V and output the identity at the end. Thus 
the decreasing permutation of size n has more than 2"" ways to be sorted i.e. more than 2" 
sorting words. 

Several articles introduce restrictions either on the rules or on the stack structure. For 
example, in his PhD-thesis West introduced a greedy model with decreasing stacks [1]. Permu- 
tations sortable with this model, called West-2-stack sortable permutations, are characterized 
and enumerated. 

For our unrestricted case called sometimes in litterature general 2-stack sorting problcTn, 
no characterization of sortable permutations and no polynomial algorithm to decide if a per- 
mutation is sortable is known. A common mistake when trying to sort a given permutation 
is to pop out the smallest element i as soon as it lies in the stacks. This operation may 
indeed move other elements if i is not the topmost element of H. The elements above it 
are then transferred into V before i can be popped out. But sometimes, it can be neces- 
sary to take some elements of a and push them onto H 01 V before this transfer. Take 



for example permutation 324617985. Trying to pop out the smallest element as soon as it 
is in the stacks leads to a dead-end. However, this permutation can be sorted using word 
P3P2>^2P4P6Pi>^ifJ'ifJ'2P7>^7>^6>^A3f^3l^4P9P8P5>^5f^5l^6l^7h^J'8>^9l^9- But we prove that this nat- 
ural idea of popping out smallest elements as soon as possible can be adapted considering 
right-to-left minima of the permutation. 

We saw that a sorting process can be described as a word on the alphabet {p, A,/i}. In 
this article, we will also describe a sorting in a different way. Take the prefix of a sorting 
word, it corresponds to move some elements from the permutation to the stacks or output 
them. At the end of the prefix some elements may be in the stacks. We can take a picture 
of the stacks and indeed, we will show that considering such pictures for all the prefixes that 
correspond to the entry of a right-to-left (RTL) minima of the permutation in H is sufficient 
to decide the sortability. Such a picture is called a stack configuration. 

Definition 1. A stack configuration c is a pair of vectors (v, w) of distinct integers such that 
the elements of v (resp. of w) corresponds to the elements of V (resp. of H) from bottom to 
top. 

A stack configuration is poppable if elements in stacks H and V can be output in increasing 
order using operations A and /u. 

Conditions for a stack configuration to be poppable have already been studied previously 
in [21 [3] and can be characterized by the following Lemma. Recall first that a permutation 
TT = 7ri7r2 . . . TTfc is a pattern of a = aia2 . . . (Tn if and only if there exist indices 1 < ii < 12 < 
. . . < ik such that ai^ai^Oi^ ■ ■ ■'^ik ^^ order isomorphic to vr. 

Lemma 2. A stack configuration c is poppable if and only if : 

• Stack H does not contain pattern 132. 

• Stack V does not contain pattern 12. 

• Stacks {V,H) does not contain pattern |2|13|. 

Moreover, there is a unique way to pop the elements out in increasing order in terms of stack 
operations. 

The first two conditions are usual pattern relation, considering elements in the stack from 
bottom to top. The third one means that there do not exist an element i in V and two 
elements j, k in H {k above j) such that j < i < k. There is a unique way to output those 
elements in increasing order as noticed in p], so we will denote by outc{I) the word that 
consists in the operations necessary to output in increasing order elements of the set of values 
/ from a stack configuration c. 

Notice that a stack configuration has no restriction upon its elements except that they 
must be different. Most of the time, a stack configuration will be associated to a permutation 
implying that the elements in the stacks are a subset of those of the permutation. In particular 
a total stack configuration of a is a stack configuration in which the elements of the stacks 
are exactly those of a. 

In this article we often use decomposition of permutations into blocks. A block B of a 
permutation a = cria2 ... (T„ is a factor iTjiTj+i . . . aj oia such that the set of values {uj, . . . ,aj} 
forms an interval. Notice that by definition of a factor, the set of indices {i, ■ ■ ■ ,j} also forms 
an interval. Given two blocks B and B' of a, we say that B < B' if and only if (Tj < aj for 



all (Tj G B, Gj G B' . A permutation a is G-decomposable if we can write \i as a = Bi . . . Bk 
such that k > 2 and for all i, Bi > -Bj+i in terms of blocks. Otherwise we say that a is 
0-indecomposable. When each Bi is 0-indecomposable, we write a = G[-Bi, . . . , -Bfc] and call 
it the 0-decomposition of a. Notice that we do not renormalize elements of Bi thus except 
-Bfc, the Bi are not permutations. Nevertheless, Bi can be seen as a permutation by decreasing 
all its elements by |-Bj+i| + ■ ■ ■ + \Bk\. 

The RTL (right-to-left) minima of a permutations are elements at such that there do not 
exist j respecting j > k and aj < a^- We denote by 0"^. the i^^ right-to-left (RTL) minima of 
a. If a has r RTL minima, then a = . . . ak^ ■ ■ ■ c^j . . . a^^ with a^-^ = 1 and k^ = n. 

Take for example permutation a = 6587413 2. The 0-decomposition of cr is a = 
[6 5 8 7, 4, 1 3 2]. Furthermore a has 2 RTL-minima which are (Jg = 1 and erg = 2. 

Definition 3. We denote a^^> = {aj \ j < ki and aj > a^-} the restriction of a to elements 
in the upper left quadrant of the i*^ right-to-left (RTL) minima a^^.. The Qi- decomposition 

of a is the Q- decomposition of a^"^' = 0[i?| , . . . , B^* ]. In the sequel Si always denote the 

number of blocks of a^'^' and B^ the j*^ block in the Qi-decomposition. 

There are two key ideas in this article. First, among all possible sorting words for a 2-stack 
sortable permutation, there always exists a sorting word respecting some condition denoted 
v. More precisely we prove that if a permutation a is sortable then there exists a sorting 
process in which the elements that lie in the stacks just before a right to left minima ki enters 
the stacks are exactly the elements of a'*'. A formal definition is given in Definition 1161 

The second idea is to encode the different sortings of a permutation respecting "P by a 
sequence of graphs Q^"^' in which each node represents a stack configuration of a block i? 
and edges gives compatibility between the configurations. The index i is taken from 1 to the 
number of right-to-left minima of the permutation. 

Section [2] study general properties on two-stack sorting and states which elements can 
move at each moment of a sorting process. Section [3] introduces the sorting graph Q^^' which 
encode all the sortings of a permutation at a given time tj and gives an algorithm to compute 
this graph iteratively for all i from 1 to the number of right-to-left minima. Last section 
focusses on complexity analysis. 

2 General results on two-stack sorting 

2.1 Basic results 

We saw that a sorting process can be described as a word on the alphabet {p, A, /i}. However 
not all words on the alphabet {p, A, p} describe sorting processes. 

Definition 4 (stack word and sorting word). Let a € {p, X, p} and w a word on the alphabet 
{p,X,p}. Then \w\a denotes the number of letters a in w. 

A stack word is a word w G {p,X,p}* such that for any prefix v ofw, \v\p > \v\x > |f|^. 

A sorting word is a stack word w such that \w\p = \w\x = \w\p. 

For any permutation a, a sorting word for a is a sorting word encoding a sorting process 
with a as input (leading to the identity of size \a\ as output). 

Intuitively, stack words are words describing some operations p, A, p starting with empty 
stacks and an arbitrarily long input and they may be some elements in the stacks at the end 



of these operations, whereas sorting words are words encoding a complete sorting process 
(stacks are empty at the beginning and at the end of the process). 

Definition 5 (subword). Let I be a set of integers. 

For any decorated word u we define u\j as the subword of u made of letters pi,Xi,iJ,i with 
i € /. For example, if u = P3tJ-5hP6P7>^6 then n|{5_6} = /isPeAe- 

We extend this definition to stack words: given a permutation a and a stack word w, wij 
is the word of {p, A, p}* obtaining from w^j by deleting indices from letters pi,Xi, pi. 

Intuitively, w\j is the subword of w made of the operations of w that act on integers of I 

Lemma 6. For any stack word (resp. sorting word) w, w\j is also a stack word (resp. sorting 
word). 

Proof. As tf is a stack word, for all i from 1 to |(t|, pi appears before Aj which itself appears 



before pi in w7j. Therefore for any prefix v of wi/, 



> \v\\ > \v\ 



If moreover i(; is a 



sorting word, let a G {p, A, p}, then for any letter Oj in w7j, pi, Aj and pi appear each exactly 



once in w7j thus \w\j\p 



\^\i\x = \W\l\^i. 



a 



Now we turn to stack configurations, beginning with linking stack words to stack config- 
urations. 

Definition 7 (Action of a stack word on a permutation). Let w be a stack word. Starting 
with a permutation a as input, the stack configuration reached after performing operations 
described by the word w is denoted Ca{w). A stack configuration c is reachable for a if there 
exists a stack word w such that c = Cij{w). In other words a stack configuration is reachable 
for a if there exists a sequence of operations p, A, p leading to this configuration with a as 
input. 

Lemma 8. If a = G[-Bi, . . . B^] then in any poppable stack configuration c reachable for a, 
elements of Bi are below elements of Bj in the stacks for all i < j (see Figure\^. 

Proof. Notice that by definition of a stack, elements of H are in increasing order from bottom 
to top for the indices. Moreover elements of V are in decreasing order from bottom to top for 
their value since from Lemma[2]they avoid pattern 12. This leads to the claimed property. D 
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Figure 2: Poppable stack configuration reachable for Q[Bi, . . . B^] 



Lemma 9. Let a be a 2-stack sortable permutation and w = uv be a sorting word for a . As- 
sume that after performing operations ofu, elements 1 ... i — 1 have been output and elements 
i . . . j are at the top of the stacks. Then there exists a sorting word w' = uu'u" for a such that 
u' consists only in moving elements i . . . j from the stacks to the output in increasing order 
without moving any other elements. 
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Proof. We claim that u' = 'U|[j...j] and u" = i'|![j...j] satisfy the desired property, where ![i ■ ■ ■ i\ 
is the set of integers [1 . . . |(t|] \ [z . . . j]. This can be checked using decorated words associated 
to w and w' and noticing that t'|[j,,,j] = ontc„(„)([z ■ ■ ■ i\) and f^[![i...j] = ■y|>j since by hypothesis 



after performing operations of u., elements 1 
at the top of the stacks. 



1 have been output and elements i . 



.J are 
D 



The stack configurations for a sorting process encode the elements that are currently in 
the stacks. But some elements are still waiting in the input and some elements have been 
output. To fully characterize a configuration, we define an extended stack configuration of a 
permutation a of size n to be a pair (c, i) where i S {l,...n + l} and c is a poppable stack 
configuration made of all elements within ui, a"2, . . . , o"j_i that are greater than a value p. The 
elements ai,...,an are waiting to be pushed and elements cjj < p,j < i have already been 
output. Notice that we don't need the configuration to be reachable. 

Definition 10. Let a be a permutation and {c,i) be an extended stack configuration of a. 
Then the extended stack configuration {c',j) of a is accessible from {c,i) if the stack config- 
uration {c',j) can be reached starting from (c, i) and performing operations p,X and /i such 
that moves /i perfomed output elements of cU {cjj . . . Un} in increasing order. 



For example, if a = 23165847 then ( LSj [6] , 7) is accessible from ( 



3t\ U, 4) by the sequence 

^,4). 



of operations fJ.2fJ'3P6P5P8^8- But (L3j L6j,5) is not accessible from 

Indeed notice that the question of whether a permutation is 2-stack sortable can be refor- 
mulated as : 
Is (i^ U, n + 1) accessible from (i^ U, 1) ? 

To solve this problem is the main goal of this article and is somehow hard, however some 
special cases are easier to deal with. The following Lemma give conditions on the involved 
configurations under which the compatibility decision problem is linear and can by solved 
by the isAccessible procedure given in Algorithm [TJ In the last sections, we show how more 
general cases can be solved using this Lemma. 

Lemma 11. Let a be a permutation of size n and {c,i), {c',j) two extended stack configura- 
tions of a with i < j. Let E be the set of elements of c and F those of d . 

• // there exists /c, ^ G {1 . . . n} such that E = {am \m <k} and F = {am \ am > ^} 

• If moreover E U F = a 

Then we can decide in linear time whether {c',j) is accessible from {c,i) using Algorithm{l\ 

Proof. We prove by case study that there is no choice between operations p, A, p at each time 
step. This is illustrated by Algorithm [T] We first prove its correctness before studying its 
complexity. 



We start with configuration curr = c. By studying specific elements of tlie current con- 
figuration curr, we prove that we can always decide which operation should be performed 
to transform curr into c' . If at any step this operation is forbidden then c' is not accessible 
from curr. Thus repeating the following process will eventually lead to decide whether c' is 
accessible from c. 

Notice that by definition, c and c' are poppable thus curr has to be poppable, hence to 
avoid the three patterns of LemmaO Let p be the next element to be output, i.e. the smallest 
element of c U {ai . . . an}- Let an (resp. cry) be the topmost element of H (resp. of V) and 
(jg be the element waiting in the input to be pushed onto H {aq may not exist and in that 
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• If ay = p then we perform fj, thanks to Lemma [H 

• Otherwise operation p, is forbidden. We have to chose between p and A. Moreover p ^V 
as V is in decreasing order from bottom to top. 

1. Suppose that an < (■■ This means that an ^ F i.e. an ^ c'. Notice that by 
definition oi p, p < an thus p ^ c' . Moreover p ^V thus p & H. U p = an then, 
by LemmalU we can pop out p. Thus we perform A. If an 7^ P, then we will prove 
that all elements x such that p < x < an form an interval at the top of the stacks. 
Those elements are all in the stacks by definition of i and p. As V is decreasing, 
the elements of [p ■ ■ ■ an] belonging to V are at the top of it. Consider now the 
position of those elements in H. 

Suppose that it is not an interval. Then it exists an element x in H such that 

X < an and there is an element y > an between x and an- But in that case, 

elements xyan form the pattern 13 2 and curr is not poppable so any movement 

is allowed here p, A or /i because we will never reach c' . 

Suppose now that the elements [p ■ ■ ■ an] form an interval in H and V. Then as 

p € H IS the smallest element, by LemmaEJ we want to pop out elements \p ■ ■ ■ an], 

hence we perform A. 

In conclusion, if an < ^ we perform A. 

2. If not, then an > ^ and thus an G c'. Once again there are different cases: 

(a) If H = then A is forbidden, thus we perform p. 

(b) If an £ H{c'), it must stay in H thus A is forbidden and we perform p. 

(c) Else an G V{c'). 

— If aq G H(c') then p is forbidden because aq would prevent an from mov- 
ing. Thus we perform A. 

— Else aq G V{c'). If aH > aq, as an G ^(c')) P is forbidden otherwise we 
cannot put aq above an in V . Thus we perform A. 

— Otherwise an-: Oq G V{c') and an < Oq. A is forbidden otherwise we cannot 
put aH above aq in V. Thus we perform p. 



We have proved that at each step of the algorithm, we know which move we have to do if we 
want to reach d . Moreover while q < j or p < i or an C ^(c')i it is impossible that curr = c' 
so we have to continue. Conversely if g > j and p > i and an ^ ^(c') then p and p and A 
are forbidden and we have to stop. Then if curr = c', c' is accessible from c, otherwise c' is 
not accessible from c. 

Finally there are at most 3n steps since at each step of the algorithm we perfom a move 
p, A or p. Moreover each step takes a constant time, therefore the algorithm runs in linear 
time. D 



Algorithm 1: isAccessible((c, i), (c', j),(t) 



Data: a a permutation and (c, i), (c' ,j) two stack configurations of a respecting 

conditions of Lemma [11] 
Result: true or false depending whether the configuration c' is accessible from c 
begin 

Put configuration c in the stacks H and V; 

p ■<— the smallest element of c U {ai . . . cr„} (next element to be output); 
q -^ i (next index of a that must enter the stacks); 
while q<j ORp<i ORaH e y(c') do 
if ay = p then 
I Perform p; p -^ p + 1; 
else 

if uh < f- then 

Perform A; 
else 

if H = OR an e H{c') then 
I Perform p; q <^ q + 1; 
else 

if aq G H{c') OR an > CTg then 
I Perform A; 
else 

Perform p; q -^ q + 1; 



Return {H, V) 



In the sequel of this article, we do not compute all possible stack configurations during 
a sorting process of a given permutation a but indeed focus on specific steps of the sorting. 
We study the possible stack configurations at each time step t, corresponding to the moment 
just before the right to left minimum ak^ is pushed onto stack H. Those configurations are 

configurations (c, fcj) accessible from (i^ U, 1). 

We will prove that we can add two different restrictions on these configurations. First, 
(c, fcj) must be a pushall stack configuration of a^^> (see below). Second {c,ki) must be an 
evolution of some configuration (c',/ci_i) between time ti-i and tj. 

Definition 12 (pushall configuration). A stack configuration is a pushall stack configuration 
of a if it is poppable, total and reachable for a. 



2.2 Prom time t, to time t. 



j+i 



Thanks to the previous decomposition into different time steps corresponding to each moment 
a right-to-left minima is pushed onto H and our previous work [3] on 2-stack pushall sortable 
permutations, we can give a polynomial algorithm deciding whether a permutation is 2-stack 
sortable. Indeed, we will prove that it is enough to consider configurations such that for each 
ti the only elements in the stacks are exactly those of ct^*-' . But a^^' is a permutation that ends 
with its smallest element such that a sorting consists in pushing all elements into the stacks 
then popping all elements out. Those possibilities are described in [3] where Proposition 4.8 
gives all possible pushall stack configurations. When a permutation is 0-indecomposable, 
Theorem 4.4 of [3] states that the number of possible pushall stack configurations is linear 
in the size of the permutation. This will ensure that our algorithm runs in polynomial time. 
Using this result, we now have the possible total stack configurations at time ti. 

The key idea for computing possible stack configurations at time ti relies on Lemma [T5l 
Informally, it is possible to decide whether a configuration at time ti can evolved into a 
specific configuration at time tj+i. Moreover, during this transition, only a few moves are 
undetermined. Indeed the largest elements won't move, the smallest one will be pushed 
accordingly to [3J and the remaining ones form a 0-indecomposable permutation that will 
allow us to exhibit a polynomial algorithm. 

First of all we denote by A^^> the common part of the permutations a^^> and a^^^^\ that 
is, A^^' = (T^*) Pi o"^*^^' = {uj \ j < ki and aj > (Jki+i}- This subpermutation A^^' intersects 
0-indecomposable blocks of a^^' and 0"^*"'""^^ Let p^^' (resp. g(*+^)) be the index such that 
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Lemma 13. For any j < min(p'-*\ (7^*"'"^)), B^ = B- 

Lemma 14. Let ai € Ay' . During a sorting process of a, elements Um such that am > cr£ 
and m < i do not move between ti and ij+i. 

Proof. Let am be an element such that m < £ and am > o'e- As a^ G A^''\ a^ > ak^_^_^ and 
j < ki, so does am > o'fci+i and m < ki. Hence both elements am, cg lie in the stacks between 
ti and tj+i (they cannot be output as ak^_^^ must be output first). Suppose that am is in H 
at time tj. As m < i, element ai is pushed after am into the stacks, thus either ai is above 
am in H or lies in V at time ti and tj+i. So, am cannot move into V, otherwise ai would be 
under it in V and V would contain a pattern 12. So, am stay in H. 

Suppose now that am is in V at time ti. As noticed previously, this element is not output 
at time tj+i. So it also lies in stack V at time t^+i, proving the lemma. D 



In the following we study conditions for 2 total pushall stack configurations c and d 
corresponding to stack configuration of o"*^*' and a^^^"^' to be accessible one from the other, 
that is, if we can move elements starting from c at time ti to obtain c' at time tj+i. 

Lemma 15. Let {c,ki) (resp. (c',A;j+i)j be a total stack configuration of a^^' (resp. a^^~^^^). 
Let vr = o", (i) . , „(»+!) then (c', /cj+i) is accessible from (c, ki) for a iff: 

1. (cf , |7r| + 1) is accessible from (c|^, tJ(-C'(*) |J -B (,)) + 1) for vr. 

2. Vj < mm(p«,g(^+i)), c (,) = c' (,). 

5. Vj > q^^'^^>,c' (j_,_i) is a reachable configuration. 



\B[ 






Proof. Suppose first that (c', /cj+i) is accessible from (c, fcj). This means that we can go from 
c to c' using operations represented by the decorated word w. These operations are stable 
that is for all /, cL is accessible from c^j. To do so, we just extract operations corresponding 
to elements of /. Indeed the decorated word wij allow to transform c into c'. This proves the 
first point of Lemma [TSl 

Let Ui G Bp . Lemma [HI ensures that elements of B^ with j < p^"^' do not move between 
ti and tj+i proving the second point of Lemma [T5l 

Finally, elements of i? for j > g(*+^) are pushed iteratively when going from c to c'. 

Those elements stay in the stacks as o"fci+i which is smaller is pushed after them. Thus they 
correspond to a pushall configuration. 

Conversely, suppose that we have the 3 different points above, we must prove that (c', A:j+i) 
is accessible from (c, ki) for a. We start by taking the stack configuration c and we will prove 
that we can obtain c' by moving elements. First of all, as c is a pushall stack configuration, 
and as elements of Bi for i > p are the smallest one and have been pushed last into the stacks 
they are at the top of the stacks (see Lemma [8|) . Thus we can pop them and output them in 
increasing order using Lemma O 

The remaining elements in the stacks don't move in the preceding operation, thus stay 

(i) 

in the same position than in c. In that configuration, elements of B\---^ are the smallest ones 
and have been pushed the latter in the stacks. Hence they lie at the top of the stacks. 

Then using point 1 of our hypothesis, we can move those elements together with pushing 
elements of B\^-_^^ \ -B (i) so that all those elements (that is elements of vr) are in the same 
position than in c'. Then, by hypothesis item 3, Vj > q^^~^^',c' ,^^1-. is a reachable configura- 



\B 



(i + l) 



tion. Thus we can push its elements into the stacks in the same relative order than in c' (see 
Lemma [8|). During these operations we ensure that elements of Bi with £ > min{p^^\ Q^^~^^')^ 
c, (i) are in the same position in our configuration than in c'. Point 2 ensures that we indeed 

obtain c'. D 

The preceding Lemma describes exactly which elements can move between tj and tj+i and 
how they move. But the hypothesis of Lemma [15] are restrictive that is configurations c and 
c' must be two total stack configurations of a^^' and (T^*"*"^^. Thus, we first prove that among 
all sortings of a 2-stack sortable permutation, there exists at least one for which the stack 
configurations at time ti contains exactly the elements of a^^' for all i. 
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Definition 16 (Properties (Pi) and (P)). Let a be a permutation and w a sorting word for 
a. w verifies (Pi) if and only if 

(i) pak.Xak.tJ'ak^ is a factor ofw. 
(a) /io-j appears before paj,_ for all aj < Cfc-. 

(Hi) All operations pai with o"£ G B^^ and j G [p^^' + l..Sj] appear before Pah+i ^^ '^■ 

where Cfc. is the i right to left minima of the permutation and o"'*-* = 0[i?| , . . . ,Bsl]. 
If a word w verifies Property (Pi) for all i then we say that w verifies Property {P). 

Lemma 17. If the sorting word encoding a sorting process of a verifies Property {Pi), then 
at time ti the elements currently in the stacks are exactly those of o"'*^ . 

Proof. By definition of time ti (just before 0"^- enters the stacks) each element in the stacks 
has an index smaller than /cj. Moreover among elements of index smaller than fcj, those of 
value greater than Ufc. cannot have been output by definition of a sorting, and those of value 
smaller than Ufc- have already been output since w satisfies item {ii) of Property {Pi). D 

Lemma 18. Let w be a sorting word for a permutation a, r be the number of RTL-minima 
of a and £ G [l--^]- If w verifies {Pi) for i G [l..i — 1] then there exists a sorting word w' for 
a that verifies {Pi) for i G [1..^]. 

Proof. Consider the sorting process of a encoded by w. The key idea is to prove that the 
smallest elements are at the top of the stacks so that we can transform the word w thanks to 
Lemma [9j 

Property (ii) for (P^) states that Haj should appear before po-^ for all aj < ak,,- Suppose 
that there still exists an element Gj with aj < a^^ in the stacks just before a^^ is pushed into 
the stacks. We prove that this element can be popped out before Ufc^ is pushed. Let aj^ be 
the smallest element still in the stacks just before Pa^- By definition, elements smaller than 
ajQ have already been output. Consider interval / = [ajQ.,aki — !]• Those elements are still 
in the stacks. If they are at the top of the stacks they can be output using Lemma [H If 
not, there exists in the stacks an element x ^ I above an element y € I. As V^ is decreasing, 
those elements are in H. Moreover x > a^^ > y- Then a^^ cannot be pushed as it will create 
a pattern 132 in H with elements x and y. Thus / is at the top of the stacks and we can 
output it before akf, is pushed onto H: using Lemma [H we build from w a sorting word w^^' 
for a satisfying {Pi) for i G [l..i — 1] and Property (ii) of {Pi). This means that w^^' can be 
decomposed as w^^' = up^^, v such that the stack configuration Ca{u) respects the following 
constraint: elements 1, . . . ,0"^^ are not in the stacks. 

So if we consider the stack configuration Ca{upa^. ), element aki is at the top of H and 

since outc^(up^ ){'^ki) = ^ak l^^k ^^ ^^^ ^^^ Lemma [9] to change the sorting word w^^' into 

a sorting word w^^' = up^^. X^^. Pa^ '^'j satisfying Property (i) for {Pi). 

Now we show considering the stack configuration c = Crj {up^^. Ao-^ p^^. ) how to transform 

the word w^"^' into a word w' = up^j^ X^^ Pa^ v^^'v^"^' with v^^> = outc{B\^i,. U • • • U Bg^). 
This will conclude the proof. 

Notice that elements of c are exactly those of a^' since the last operations performed are 
/9(jfe Ao-j. Pak ^iid elements are pushed in the stacks in increasing order of indices and output 
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in increasing order of values. Thus outc{B \^^i. U • • • U Bs/) = out{Bsi ) ■ ■ ■ out{B\^^s ) (see 
Lemma [8]). We show by induction on j from S£ to p'*^ + 1 that we can build a sorting word 
for a of the form up^^ X^^ p.^^. v^^'^'v^'^'^' with v^^'^' = out{Bs/) ■ ■ ■ out{B- ). For j = si that 
is a word in which elements of block Bg^ are output immediately after 0"^^ has been output. 
By definition of sg and because elements of c are exactly those of o"' \ all elements of Bs^ 
lie in the stacks in configuration c, are the smallest elements in this configuration and lie at 
the top of the stacks in configuration c (see Lemma [8]). Hence, using Lemma |9], there exist a 
sorting word w^^> for a such that w^^'' = up^j^, Aq-j. p-a^, out{BsJv". Repeating this operation 

for all blocks Bj with j from S£ — 1 to p^'^' + 1, we have Property (iii). D 

Notice that Property (Pq) is an empty property satisfied by any sorting word. Using 
recursively Lemma [18] we can transform any sorting word into a sorting word satisfying 
Property (P), leading with Lemma [17] to the following theorem: 

Theorem 19. If a is 2-stack sortable then there exists a sorting word of a respecting Property 
(P) . In particular, in the sorting process that this word encodes, the elements currently in the 
stacks at time ti are exactly those of a^^' . 

Theorem 1 1 91 ensures that if a permutation is sortable then there exists a sorting in which at 
each time step ti, elements in the stacks are exactly those of 0"'*^ Thus stack configurations at 
time ti and tj+i satisfy hypothesis of Lemma [T5] and we can apply it to decide if a permutation 
if 2-stack sortable. 

3 An iterative algorithm 

3.1 A fisrt naive algorithm 

From Theorem [19] a permutation a is 2-stack sortable if and only if it admits a sorting 
process satisfying Property (P). The main idea is to compute the set of sorting processes of 
a satisfying Property (P) and decide whether a is 2-stack sortable by testing its emptiness. 

Verifying (P) means verifying (Pj) for all j from 1 to r, r being the number of right-to-left 
minima (whose indices are denoted kj). The algorithm proceeds in r steps: for all i from 1 
to r we iteratively compute the sorting processes of cr<fcj verifying (P^) for all I from 1 to i. 
As a<kr = '^) the last step gives sorting processes of a satisfying Property (P). 

By "compute the sorting processes of CF<_ki^ we mean compute the stack configuration 
just before Ufc. enters the stacks in such a sorting process. Note that this is also the stack 
configuration just after at has been output since p^^ X^^ Paj, is a factor of any word verifying 
(P). 

Definition 20. We call Pi-stack configuration of a a stack configuration Ca{w) for which 
there exists u such that the first letter of u is pa,.. and wu is a sorting word of a<:k. verifying 
(P) for cr<ki (that is, verifying (P^) for all i from 1 to i). 

Lemma 21. For any i from 1 to r, cr<j^. is 2-stack sortable if and only if the set of Pi -stack 
configurations of a is nonempty. In particular, a is 2-stack sortable if and only if the set of 
Pr-stack configurations of a is nonempty. 

Proof This is a direct consequence of Definition [20] and Theorem [19] D 
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Lemma 22. Any Pi-stack configuration of a is a pushall stack configuration of a^^' accessible 
from some Pi^i-stack configurations of a. 

Proof. By definition of (P), each Pj-stack configurations of a is accessible from some Pi^i- 
stack configurations of a (take the prefix of w that ends just before po-^. )• Moreover it is a 

pushah stack configuration of cr^*' from Lemma [T71 D 

As explained above, the algorithm proceeds in r steps such that after step i we know 
every Pj-stack configuration of a and we want to compute the Pj+i-stack configurations of 
a at step i + 1. As configurations for i + 1 are a subset of pushall stack configurations of 
a^^~^^' , a possible algorithm is to take every pair of configurations {c,c') with c being a Pi- 
stack configuration of a (computed at step i) and c' be any pushall stack configuration of 
^(«+i) (given by Algorithm 5 of [3]). Then we can use Algorithm [1] to decide whether c' is 
accessible from c for a. This leads to the following algorithm deciding whether a permutation 
a is 2-stack sortable: 
Algorithm 2: isSortableNaive 
Data: a a permutation 

Result: true or false depending whether a is 2-stack sortable 
begin 

E, F two empty sets; 
E^ PushallConfigs(fT(i)); 
for i from 2 to r do 
F ^ 0; 
for c in E do 

for c' in PushallConfigsfa^"^' ) do 

if isAccessible({c, hi), (c', ki^i),a) then 
\_ F^FUc'; 

E^F; 

if E is empty then 
I return false; 
else 

return true; 

Notice that at step i, the set E computed contains all Pj-stack configurations of a but 
may contain some other configurations. However since each configuration of i? is a pushall 
configuration of a^"^' and is accessible for a from some pushall configurations of a^'^~^\ each 
configuration of E indeed corresponds to some sorting procedure of i7<fc. , proving the correct- 
ness of Algorithm [2l 

But this algorithm is not polynomial. Indeed the number of Pj-stack configurations of a 
is possibly exponential. However this set can be described by a polynomial representation as 
a graph Q^^> and we can adapt Algorithm [2] to obtain a polynomial algorithm. In this adapted 
algorithm, the set E computed at step i is exactly the set of Pj-stack configurations of a. 
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3.2 Towards the sorting graph 

We now explain how to adapt Algorithm [2] to obtain a polynomial algorithm. Instead of 
computing all Pj-stack configurations of u (which are pushall stack configurations of o"'-*)), we 



compute the restriction of such configurations to blocks 5- of the Q-decomposition of u*-*'. 



J 

By Lemma [HI those configurations are stacked one upon the others. The stack configurations 
of any block B- are labeled with an integer which is assigned when the configuration is 
computed. Those pairs configurations / integer will be the vertices of the graph ty(*) which we 
call a sorting graph, the edges of which representing the configurations that can be stacked 
one upon the other. Vertices of the graph Q^^' are partitioned into levels corresponding to 
blocks By . To ensure the polynomiality of the representation, we will prove that a given 
integer label could only appear once per level of the graph Q'^^K As those numbers are 
assigned to configurations when they are created, each integer corresponding to a pushall 
stack configuration, from [3] there exists only a polynomial number of disctincts integers thus 
of vertices. This will be explained in details in the next section. The integer indeed can be 
seen as the memory of the configuration that encodes its history since it has been created: two 
configurations which have the same label come from the same initial pushall configuration. 

More precisely a sorting graph Q'^^' for a permutation a of size n and an index i verifies 
the following properties: 

• Vertices of Q^^> are partitioned into Sj subsets V with j £ [I . . . Si]. 

• For any j E [1 . . . Si], IVJ^ | < 9n + 2. 

• Each vertex v € Q^"^' is a pair (c, ^ ) with c a stack configuration and ^ an index 
called configuration index. 

• All configuration indices are distinct inside a graph level V- 

• (C) ® ) € V ^ c is a pushall stack configuration of B- accessible for a. 

• There are edges only between adjacent blocks V , V^^-^. 



• 



(i) (i) ("\ 

Paths between vertices of V^ and Vsi corresponds to stack configurations of a^"^' . More 
precisely such paths are in one-to-one correspondence with Pj-stack configurations of 
a (that is, stack configurations corresponding to a sorting of (T<fc. respecting (P) just 
before 0"^- is pushed onto H). 



V. 



I ') (i) (i) 

For any vertex v of t/"', there is a path between vertices of V^ and Vs/ going through 



Though the definition of sorting graph is complex, its use will be quite understandable 
and easy. Look for example at the permutation a = 4321. There is only one right to left 
minimum which is 1. Compute all possible stack configurations just after 1 enters H. At this 
time, all elements are in the stacks since the first element which must be output is 1. More 
formally, we are looking at the pushall stack configurations of a with 1 in H. 

There are 8 different such configurations which are: 
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^ 



^1 



w 



m 



The 0-decomposition oi a is a = 0[4, 3, 2, 1]. We build a graph with 4 levels, each level 
corresponding to pushall stack configurations of a block. 




Stack configurations of i?i = 4 



Stack configurations of i?2 = 3 



Stack configurations of B^ 



Stack configurations of -B4 = 1 



Figure 3: Graph encoding pushall stack configurations of a = 4321. 



Then the 8 configurations are found taking each of the 8 different paths going from any 
configuration of Bi to configuration \_^\l\ of S4. In Figure El the thick path gives the stack 

2 1 

configuration [^ [3j by stacking the selected configuration of ^4 above the configuration of 

B3 and so on. 

But in the last level B4 we only consider configuration 1^ LlJ so this level is useless. The 
sorting graph Q^^' for a = 4321 encodes pushall stack configurations of a^^' = 432, corre- 
sponding to stack configurations just before 1 enters H (and not after as above). 

There are 8 different such configurations which are: 



m 



m 



4jy y^ 



m 



As the G-decomposition of a^^' is a^^' = 0[4, 3, 2], the sorting graph Q^^> has 3 levels. 



15 




Stack configurations of i^i = 4 



Stack configurations of i?2 = 3 



Stack configurations of S3 = 2 



Figure 4: Sorting graph ^'-^^ of o" = 4321. 



Then the 8 configurations are found taking each of the 8 different paths going from any 
configuration of Bi to any configuration of ^3. In Figure IH the thick path gives the stack 



[3j by stacking the selected configuration of -63 above the configuration of 



configuration 
B2 and so on. 

We transform Algorithm [2] to a polynomial algorithm by computing at step i not all Pi- 
stack configurations of a, but instead the sorting graph Q^'^' encoding them. The graph Q^'^' 
is computed iteratively from the graph Q^'^~^' for any i from 2 to r. The way G^^' is computed 
from Q^^~^' depends on the relative values of p^*' and g'*"*"^'. By definition of a sorting graph 
given piMl if at any step G^^' is empty, it means that (T<fc. is not sortable (from Theorem I19p 
and so is a thus the algorithm returns false. This is summarized in Algorithm [31 

Algorithm 3: isSortable 
Data: a a permutation 

Result: true or false depending whether a is 2-stack sortable 
begin 

Q ^ ComputeGl; 
for i from 2 to r do 
if p(i) = g('+i) then 
I G ^ iteratepEqualsq{G) or return false 
else 

if pW < g(^+i) then 

I G ^ iteratepLessThanq{G) or return false 
else 
\_ G ^ iteratepGreaterThanq{G) or return false 

return true 
In the next subsections we describe the subprocedures used in our main algorithm isSortable{a). 



3.3 First step: ^(^) 

In this subsection, we show how to compute the Pi-stack configurations of a, that is, the 
stack configurations corresponding to time ti for sorting words of cr^k^ that respect (P) for 

From Lemma [22I such a stack configuration is a pushall stack configuration of a^^' . Con- 
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versely since ak^ = 1, cr^^' = (y<kx and each sorting word oia<,kx respects (Pi) for cr<fc^. Thus 
the set of Pi-stack configurations of a is the set of pushall stack configurations of a^^' . 

By Proposition 4.7 of [3j, these stack configurations are described by giving the set of 
stack configurations for each block of the G-decomposition of a^^> . More precisely, with cr'^' 



?(i) 



?(i) 



\B\ , . • • ) Bsi] there is a bijection fiom pushaUConfigs{B\ 



(1)^ 



(1)^ 



xpushallConfigs{Bsi 



onto pushallConfigs{a^^') by stacking configurations one upon the other (as in LemmalH]). As 
a consequence, from Lemma [21] (T<fc^ is not sortable if and only if a set pushallConfigs{B- ) 
is empty. 

Moreover it will be useful to label the configurations computed so that we attach a distinct 
integer to each stack configuration when computed. 

At this point, we have encoded all configurations corresponding to words respecting P up 
to the factor piXifii. 

The obtained graph is Q^^' . This step is summarized in Algorithm |4l 
Algorithm 4: Computed 

Data: a a permutation, num. a global integer variable 

Result: false if a<ki is not sortable, the sorting graph Q^^' otherwise. 

begin 
E = 0; 

Compute a^^> and its 0-decomposition Q[B[ , 

for j from 1 to s\ do 

Vj'^ ^ 0; 



,B. 



(1)1. 

Sl J) 



S = pushallConfigs{B- 
if 5 = then 

return false; 
else 

for s € 5 do 
(1) 



(1)^ 



V 



vP[]{{s, 



] 



1 
num. ^ num + 1; 

if J > 1 then 




)}; 



I E = E[j{is,s'),s ev}'\s' evl'_\} 



return e(^) = ( U V} ,E) 



je[i..s[ 



{i)i 



3.4 From step i to step i + 1 

After step i we know the graph Q^"^' encoding every Pj-stack configuration of a and we want 
to compute the graph ^'*+^) encoding Pj+i-stack configurations of a at step i + 1. From 
Lemma [22] we have to check the accessibility of pushall stack configuration of a^^^^' from 
Pj-stack configurations of a. We want to avoid to check every pair of configurations (c, c') 
with c being a Pj-stack configuration and c' be a pushall stack configuration of a^^'^^^ because 
the number of such pair of configurations is possibly exponential. Thus our algorithm focuses 
not on stack configurations of some a^^^ but on sets of stack configurations of blocks P- , 
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making use of Lemma [T5l Using Lemma [22l Lemma [15] can be rephrased as: 

Lemma 23. Let d he a total stack configuration of a^"^^^' , p = p^^' and q = 5*-*+^'. Then c' 
is a Pi+i-stack configuration of a if and only if: 

• For any j < q, c' (^^i) is a pushall stack configuration of a. a+i), and 

• There exists a Pi- stack configuration c of a such that : 



,{») 



I iTiinfp.g) ' 



(j) is accessible from c. (i+i) 



inin(p,q) t^ \ P '^J h 



Recall that a Pj-stack configuration of a is encoded by a path in the sorting graph ^(*', 
corresponding to the 0-decomposition of the permutation 0"^*^ into blocks B^ . The last point 
of Lemma [23l ensures that the first levels (1 to min{p^^' , q^'^'^^') — 1) are the same in ^(*+^) 
than in Q^^' . The first point of Lemma [23] ensures that the last levels (> q^^~^^') of ^(*+^) form 
a complete graph whose vertices are all pushall stack configurations of corresponding blocks. 
So the only unknown levels for ^(*+^) are those between min(p'*\g(*"'"^)) and q^^^^' and we 
can compute them by testing accessibility. 

There are differents cases depending on the relative values of p^*^ and q^'^'^^' . To lighten 
the notations in the following, we sometimes write p (resp. q) instead oi p^''> (resp. q^^~^^>). 



3.4.1 Case p(*) = g(*+i) 



Upii) = g(i+i) then B\J^l n A(^) = B% n A^'^ (see Figure [5]) 



(i) 



□ 




B 




'l; 






















0"fc.- 


o-t, 






□ 
□ 





D 



□ 



^^,+1 B^;; 
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'^k. 



D 



□ 



□ I 






n 



□ 



D 



□ 

• 



ffci 



Figure 5: Block decomposition of a^^' and of a^^^^' when p^^' = g(*+^) 



We have the sorting graph Q^'^' encoding all Pj-stack configurations of a and we want to 
compute the sorting graph ^(*+i) encoding all Pj+i-stack configurations of a assuming that 

In this case, from Lemma [23] we only have to check accessiblity of pushall configurations 



P 



of Bq from configurations of Bp belonging to level p of Q^'^K Indeed from the definition 
of a sorting graph given pll4[ for any vertex v of ^'*^ there is a path between vertices of V^ 



(i) 



r{i) 



and Vsi going through v, and such a path corresponds to a Pj-stack configuations of a. Thus 
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for any configurations x of Bp belonging to a vertex v of level p of Q^'''\ there is at least one 



Pj-stack configurations c of o" such that c^ {i) 

I ^p 



X, and c.^{i) (i) is encoded by a 

I 1 min(p,q) — 1 



path from v to level p of Q^"^' (which go through each level < p). 

If there is no pushall configuration of Bq accessible from some configurations of Bp 
belonging to level p of G^^, or if o"*^*^^' has no pushall configuration, then a has no Pj+i-stack 
configuration and cr<ki+i is not sortable (from Lemma [2T]l . 

This leads to the following algorithm: 

Algorithm 5: iteratepEqualsq{Q 



(i)^ 



Data: a a permutation and G^^' the sorting graph at step i 

Result: false if 0"<fc.^^ is not sortable, the sorting graph (j(*+^) otherwise. 

begin 

G an empty sorting graph with Sj+i levels; 

G' ^ ComputeGl{a^'^^^') (pushall sorting graph of a^"^^^) or return false; 

Copy levels q + 1 . . . Sj+i of G' into the same levels of G'-, 

for (c, ^) in level p of G^^^' do 

Ti the subgraph of G^^' induced by (c, ©) in levels < p; 
for (c',@) in level q of G' do 

sAccessible{c,c',a. (i). , „{»+!)) then 

Add (c',^p) in level q oi G (if not already done); 
Merge Ti in levels < g of ^ with (c',^p) as origin; 



if level q of G is empty then 
I return false; 

for (c',^) in level q of G do 

I Add all edges from (c', ^) to each vertex of level g + 1 of ^; 
_ return G 



3.4.2 Case p(*) < g(*+i) 

lipii) < qii+^) then S^!+V) n ^« ^ B^!|) n yl» (see Figure ED- 




W. — ^n 



o(~i) 



□ 



<^ki 
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□ 
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^"^^ — ^n 
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'^fe,+i 
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Figure 6: Block decomposition of a^^' and of cr'-*^^-' when p*-*-* < g(*+i) 
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Again, Lemma [23] ensures that the first p — 1 levels of (y(*+^) come from those of G^^' and 
the levels > q are all pushall stack configurations of the blocks B^ ' of a^"^^^' . The difficult 
part is from level "p to level q. As in the preceding case, by Lemma [23l we have to select 
among pushall stack configurations of blocks |?,p + 1, . . . , g of a^'^^^' those accessible from a 
configuration of B^ that appears at level p in (/(*). We can restrict the accessibility test 
from configurations of Bp appearing in graph C/(*) to pushall stack configurations of Bq . 
Indeed, Lemma O ensures that elements of blocks i? for j from p to g — 1 are in the 

same stack at time tj and at time tj+i. Thus configurations of i? for j from p to g — 1 

are restrictions of configurations oi Bp . We keep the same label in the vertex to encode that 
those configurations of Bp , -Bp^i ' • • • > -^o-i come from the same configuration of Bp and 
we build edges between vertices of -By^i ^^id -^i* that come from the same configuration of 
Bp . It is because of this case p = q that we have to label configurations in our sorting graph. 
Indeed two different stack configurations ci and C2 of Bp may have the same restriction to 
some block B^^ but not be compatible with the same configurations, thus we want the 
corresponding vertices of level j of ^'*+^) to be distinct, that's why we use labels. 

More precisely we have the following algorithm. 

Algorithm 6: iteratepLessThanq[Q^^>) 

Data: a a permutation and Q^^' the sorting graph at step i 

Result: false if o"<fc.^^ is not sortable, the sorting graph Q'^^'^^' otherwise. 

begin 

G an empty sorting graph with Sj+i levels; 

Q' ^ ComputeGl{a^^^^') (pushall sorting graph of a^"^^^') or return false; 

Copy levels q + 1, . . . , Sj+i of Q' into the same levels of G; 

for (c, ©) in level p of Q^^' do 

Ti the subgraph of Q^^' induced by (c, ®) in levels < p; 
for (c',^) in level q of Q' do 

if isAccessible{c,c',a.(i). .„(i+i)) then 

Add (c',^p) in level q oi Q (if not already done); 
for j from, q — 1 downto p do 



Add (C|„(i+i),®) in level j of Q; 

I 7 



Add an edge between (c 



1^(^+1) 



,0) and (C| (i+i),©) in Q. 



Merge Ti in levels < p oi G with (c, „(i+i) , Qi) as origin; 



if level q of Q is empty then 
return false; 

for (c',^p) in level q of Q do 
I Add all edges from (c',^) to each vertex of level q + 1 oi Q; 

return Q; 



Note that in Algorithm [6l before calling isAccessible{c,c' ,a. 



B^'>U4'+''- 



we extend con- 



figuration c' to D^"^' U Bq^ ' by assigning the same stack than in c to points of D^'^' \ Bq ' . 
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This is justified by Lemma [T^ 
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3.4.3 Case p^'^ > g(*+i) 

lipii) > qii+i) then B% n yl» £ B^\t^l n yl» (see Figure ED- 
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Figure 7: Block decomposition of cr'*-' and of cr'*"'"^-' when p'*-* > g(*+^) 



Algorithm 7: iteratepGreaterThanq{Q^'^') 



Data: a a permutation and Q^^^ the sorting graph at step i 

Result: false if <7<ki+i is not sortable, the sorting graph Q^^~^^> otherwise 

begin 

Q an empty sorting graph with Sj+i levels; 

Q' ^ ComputeGl{a^^^^') (pushall sorting graph of o"'*^^^) or return false 
Copy levels q + 1, . . . , Sj+i of Q' into the same levels of Q; 
for (c, @) in level p of Q^^' do 
for (c',^M) in level q of Q' do 

if isAccessible(c,c'a,„(i), i„{i+i)) then 

\t!p \JBq 

if there is a path (c, O) o (c' (,, , ^) o . . . o (c' (,, , 
Add (c',^) in level q of Q (if not already done); 



') in C/(*) then 



v. the subgraph of G^^' induced by (c' {,),^™) in levels < q; 

\Bq ^1^ 

Merge Ti in levels < g of ^ with (c',^) as origin; 



if level q of Q is empty then 
return false; 

for (c',^) in level q of Q do 
I Add all edges from (c',^p) to each vertex of level g + 1 of ^; 

return Q; 

(i) ■ 

This case is very similar to the preceding one except that Bp is not cut into pieces 
but glued together with preceding blocks. As a consequence, when testing accessibility of a 
configuration of Bg , we should consider every corresponding configuration in Q^"^' , that is 
every configuration obtained by stacking configurations at level q,q + 1, ... ,p in Q^^K Unfor- 
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tunately this may give an exponential number of configurations, but noticing that by Lemma 

9+1 •• • ^P- 



[TH elements of blocks Bq,B^V^ . . . -Bp_^ are exactly in the same stack at time tj and at time 



tj+i, it is sufficient to check the accessibility of a pushall configuration c' of Bq from a 
configuration c oi Bp and verify afterwards whether the configuration c has ancestors in Q^^' 
that match exactly the configuration c'. This leads to the Algorithm [71 

Note that in Algorithm [71 before calling isAccessible{c,c',a.{i). .^(i+i)) we extend con- 

figuration c to Z?'*^ |J Bp by assigning the same stack than in c' to points of D^^' \ Bp . This 
is justified by Lemma [TH 

Now that we have described all steps of our algorithm, we turn to the study of its com- 
plexity. 

4 Complexity Analysis 

In this section we study the complexity of our main algorithm: isSortable(a) (Algorithm [3]) . 
The key idea for the complexity study relies on a bound of the size of each graph Q^^' , as 
described in the following lemma. 

Lemma 24. For any i € [^■■t], the maximal number of vertices in a level of Q^^> is 9n + 2 
where n is the size of the input permutation. 

Proof. From Theorem 4.4 of p], the maximal number of pushall stack configurations of a 
0-indecomposable permutation vr is 9|7r| + 2. 

By definition of Q^^' , the vertices of a level correspond to pushall stack configurations of 
a given block of the ©i-decomposition of the input permutation a. Thus the cardinality of 
a level is bounded by 9A; + 2 where k is the size of the corresponding block. As A; < n, the 
result holds for i = 1 . 

Suppose now that the result is true for a given Q^'^' , we show that it is then true for ^(*+^). 
The graph ^(*+^) is build from Q^^' using Algorithm [H [H or [71 In each case for a level j of 
^(«+i) we have: 

If j > g(*+^) then vertices of the level j of (/(*+^) are the pushall stack configurations 
corresponding to the block B^^ of the ©j+i-decomposition of a. Thus Theorem 4.4 of [3] 
ensures that the cardinality of level j is bounded by 9n + 2. 

If j = g(*+^) then vertices of the level j of ^(*+^) are a subset of the pushall stack config- 
urations corresponding to the block B^^ of the 0j+i-decomposition of a. Again Theorem 
4.4 of p] ensures that the cardinality of level j is bounded by 9n + 2. 

If j < p^^' then vertices of the level j of Q^^~^^) are a subset of vertices of the level j of Q^^' . 
By induction hypothesis the cardinality of level j is bounded by 9n + 2. 

If p(^J < j < q^^~^^' then vertices of the level j of ^(*+^) are restrictions of a subset of 
vertices of the level j of Q^"^' . By induction hypothesis the cardinality of level j is bounded 
by 9n + 2, concluding the proof. D 

Lemma 25. For any i € [l..r], the number of vertices of Q^''"' is 0{n'^) and the number of 
edges of Q^''' is 0{n^), where n is the size of the input permutation. 

Proof. The result follows from Lemma [2H as there are at most n levels and there are edges 
only between consecutives levels. D 
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Theorem 26. Given a permutation a, ^/^orit/imO isSortable((T) decides whether a is sortable 
with two stacks in series in polynomial time w.r.t. \a\. 

Proof. Algorithm [3] involves four other subroutines: Computed (Algorithm S]), iterate- 
pEqualsq (Algorithm [5]) , iteratepLessThanq (Algorithm [6]) and iteratep Greater Thanq (Algo- 
rithm [7]). 

Each for-loop in these algorithms is executed at most a linear number of time by Lemma 

m 

Moreover each included operation is polynomial by Lemmas [25] and [TT] D 

Notice that a more precise analysis of complexity leads to an overall complexity of O(n^). 
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